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Abstract 

We suggest a new adiabatic approach for description of three charged 
particles in the continuum. This approach is based on the Coulomb- 
Fourier transformation (CFT) of three body Hamiltonian, which al- 
lows to develop a scheme, alternative to Born-Oppenheimer one. The 
approach appears as an expansion of the kernels of corresponding in- 
tegral transformations in terms of small mass-ratio parameter. To be 
specific, the results are presented for the system ppe in the continuum. 
The wave function of a such system is compared with that one which 
is used for estimation of the rate for triple reaction p + p + e — > d + u, 
which take place as a step of pp-cycle in the center of the Sun. The 
problem of microscopic screening for this particular reaction is dis- 
cussed. 

1 Introduction. 

An accurate treatment of three charged particles in the continuum 
at low relative energies represents till now a very difficult problem 
which is actual in many areas of physics. Indeed, only some special 
cases with specific properties of the system has been considered in 
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the literature. In this respect one should mention papers [1, 2] where 
the asymptotic solution for the three-body wave function has been 
obtained for configurations when all interparticle distances are much 
larger of the characteristic size of the system. An alternative limiting 
case considered in [3] corresponds to configurations where one Jacobi 
coordinate is much larger than the other one. The near threshold 
breakup of hydrogen by proton (or electron) studied in [4] is another 
example of approximative solution of three body Schrodinger equation 
obtained so far for three charged particles in the continuum. 

The purpose of this paper is to develop a new adiabatic expansion 
for three-body Hamiltonian for the system consisting of one light and 
two heavy particles. The asymptotic behavior of the wave-function 
with respect to the coordinate of light particle will be presented in 
the framework of this new adiabatic expansion. We use the Coulomb- 
Fourier transform formalism proposed in [5] to make a unitary trans- 
formation of the Hamiltonian, which leads to a convenient represen- 
tation, where, for example, one long-range interactions is eliminated. 
The known explicit form of eigenfunctions allows us to construct the 
useful integral representation of interaction potentials for transformed 
Hamiltonian. The important feature of this representation for poten- 
tials is the appearance of universal integral with the integrand contain- 
ing the exponential factor, proportional to square root of mass ratio 
r ~ ^Jm e /m p of light (electron) and heavy (proton) particles. The 
natural power series in r of this integral then generate power series 
representation of the transformed Hamiltonian what is the basis for 
our new variant of adiabatic representation of the problem. As the 
first stage, in this paper we study the solution of the problem taking 
into account terms up to 0(r 2 ) order. 

The paper is organized as following. Section 2 plays the central 
role and includes formulation of the problem, description of necessary 
portions of Coulomb-Fourier transformation. Here we give as exact 
formulas for CF transformed Hamiltonian for ppe system as well as 
derive the expansion of matrix elements of the Hamiltonian and wave 
function as power series in r 2 . Some technical details concerning eval- 
uation of singular integrals from the main tex is placed in Appendix. 
Third section contains application of approximation for three-body 
wave function obtained in Section 2 to description of some reactions 
of pp cycle. Short conclusion summarizes the paper. 

In the paper we use units such that fi = c = 1, for the unit electric 
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charge the symbol e is used. Three-dimensional vectors are denoted by 
x, y, k, p... and its modules and respective unite vectors by x, y, k,p... 
and by x, y, k, p... . Sometimes we combine pairs of three-dimensional 
vectors in six-dimensional ones as X = {x, y}, P = {k, p}. The 
Hilbert space of functions depending on vectors X which in our paper 
play the role of configuration space vectors will be denoted as H and 
the Hilbert space H will be associated with function depending on 
momentum variables P. 

2 Adiabatic expansion for three-body 
Hamiltonian and solution 

The ppe system where p is proton and e is electron with masses m p 
and m e (m p S> m e ) respectively is considered. Number 1 we assign to 
electron whereas 2 and 3 to protons. The Hamiltonian of the system in 
the center of mass frame using mass-renormalized Jacobi coordinates 
can be written in the form 

H = -A Xl - A yi + v s (x 1 ) + ni/xi + n 2 /x 2 + n 3 /x 3 . 

Here V s (x\) is a short-range potential describing strong pp interaction. 
Mass-renormalized charge factors rii are defined by the formulas rii = 
ejeky/2[ijk, where e\ = — e, e 2 = = e are electron and protons 
charges and iiij stands for reduced mass of a pair of particles ij, i.e. 
fiij = mimj/(rrii + rrij). Introducing proton and electron masses into 
this formula we get ^23 = m p /2, ^31 = m e m p /(m e + m p ), [112 = //31. 

Before proceed further, let us make three clarifying comments. 
First, throughout the paper we systematically use a convention that 
indices of any pair of particles ij are considered as a part of triad ijk 
which itself is a cyclic permutation of 1, 2, 3. Second, let us recall that 
mass-renormalized Jacobi set Xj, y^ is defined in such a way that the 
vector Xj up to the factor y / 2fijf c is proportional to the relative position 
vector of particles j and k and the vector y« is the position vector of 
particle i with respect to the center of mass of corresponding two- 
body sybsystem. There are three possible sets x«, yj, i = 1,2,3 and 
different sets are related to each other by kinematic rotation relations 

Yi = -Sij^i + Cijyj (1) 
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with coefficients being denned in terms of particie masses by the for- 
mulas 

• ' (2) 



3 \/ {mi + m k ){mj+m k y 

Sij = eijk^Jl — cfj, where eijfc is fully antisymmetric tensor normal- 
ized as £123 = 1. Third, in all mass factors we keep general sit- 
uation, i.e. not neglecting m e with respect to m p , making some- 
times simplifications for illustrative purposes, as for example ^31 = 
m e m p / (m e + m p ) = m e (l + 0(m e /m p )). 

Now let us turn to the solution of the Schrddinger equation 

= E^> (3) 

for three particles in the continuum {E > 0). To this end we will 
construct a special representation. As the basis of this representation 
we take the eigenfunctions \I/ c o(X, P), i.e. H c o^ c q = P 2 ^ c o, of the 
operator 

H c0 = -A X1 - A yi + ni/xi 

with repulsive (n\ > 0) Coulomb potential. It is obvious that this 
eigenfunctions have the form 

^o(X,P)=^ 1 (x 1 X(y 1 ). 
Here V'pi(yi) = ( 27r )3/2 e^ Pl,yi ^ is normalized plane wave and 

^(xi) = T2^e i < kl ' xl >e-^/ 2 r(l + i 7 i)*(-i7i,l,ifciei) (4) 

is the normalized Coulomb wave function. The standard notations 
for Sommerfeled parameter 71 = n\/2k\, parabolic coordinate £1 = 
xi — (xi,ki), Gamma function T(z) and Confluent Hyper geometric 
function 5>(a, c, z) have been used. 

The representation which we call Coulomb-Fourier (CF) one is 
generated by the transform [5] 

G(X) = T C G = J dP * c0 (X, P)G(P). (5) 

The integral operator J- c transforms the Hilbert space 7i of functions 
depending on momentum variables P into Hilbert space TC of functions 
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depending on coordinates. Note that in the limiting case n\ = 
the T c operator is reduced to the standard inverse Fourier transform 
which connects conventional momentum space and configuration space 
representations. The Hilbert space adjoint T\ acts from Ti to 7i and 
for the pair T c and T\ the unitarity properties hold 

J~\j~ 2 = Iy^, J~ C J~\ = I'M, 

which are just the operator form of the orthogonality and completeness 
of eigenfunctions of the Hamiltonian H c q. 

In the Schrodinger equation in CF representation described above 

H* = EV, (6) 

the CF transformed Hamiltonian H = T\RT C appears now as integral 
operator with the kernel (matrix elements) 

^(P,P') = «^ 1 W^ i ^ i )= (7) 

(k?+p?)«5(k 1 -k' 1 )«5(p 1 -pi)+^(k 1 ,k' 1 )<5(p 1 -pi)+W 2 (P, P')+W 3 (P, P') 

operating on CF-transformed wave function ^(P). Here the first term 
corresponds to kinetic energy operator Hq = kf + pf , in the second 
term v s stands for the CF-transformed short-range pp interaction po- 
tential 



^(ki,ki) = (Y> kl HY> k; ) = J <&iC( x i) u -( x i)^ki( x i) 

where * means the complex conjugation, and Wj are Coulomb poten- 
tials rij/xj in CF representation. Let us notice, that the contribution 
from n\jx\ potential has been eliminated by CF transform. The func- 
tions Wj(P,P'), j = 2,3 have the following structure 

w^(p,p') = i^iI-^^Hp-pOK^p.p'). 

where 

is the familiar Fourier transform of Coulomb potential rij/xj and the 
functions £j(P,P'), j = 2,3 are given by the integrals 

£j(P,P')= lim /dx 1 e^< Xl ' p - p '>- A l Xl IV'k* 1 (xi)V'k'( x i)- ( 8 ) 

A — ►+() J 1 
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The parameters Tj, j = 2,3, have the kinematical origin and are rep- 
resented in terms of kinematic rotation matrix elements (1) as 

Tj = Cjl/Sjl. 

Noting that T3 = — r 2 which is the consequence of equality of heavy 
particles (protons) masses and using definitions (2) and the fact that 
m e <C m p we get 

r 2 = y / m e /2m p (l + 0(m e /m p )), 

what shows that Tj are small. This allows us to expand the exponential 
factor in the integrand of (8) and obtain in general case the expression 

£,(P,P') = 5(ki-k / 1 )+ (9) 



2! 

Here L^(P,P') are integrals 



n LW (P, P') + ^U (2) (P, P') + i^L( 3 ) (P, P') 



L«(P,P') = lim / dx ie - A l Xl IC(xi)(x 1 ,p 1 -p / 1 )V^ 1 (xi) (10) 

which are independent on j. This last fact and the property T3 = — r 2 
leads us to the following expansion for the sum of the CF transformed 
Coulomb potentials W2 + W3, which contains only even power terms 

W 2 (P,P')+W 3 (P,P')=V C eff ( Pl ,p' 1 )x (11) 

{<5( kl - k' x ) - ^L( 2 )(P,P') + ^W(P,P') - ...} , 

where we have introduced parameter r = |t2 | - The quantity £>gjj(pi, p^) 
stands for the Coulomb potential corresponding to interaction between 
electron and effective particle with the charge 2e and the mass 2m p 
and has the form 

%fM^ = ^A («) 



with n e ff = -2e 2 \/2m e ~ -2e 2 ^2 



m e 2m p 
m e +2m p ' 
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Inserting formula (11) into (7) we arrive at the representation of 
the CF-transformed Hamiltonian H which plays a central role in the 
solution of the problem 

H = H + V S + V e c ff + t 2 W. (13) 

The kernels of operators involved in (13) read 

flo(p,p') = (ki + Pi) <y(p-p'), 

£(P,P') =fi s (ki,ki)<5( P i-p , i), 
K c // (P,P')=t) e c // (pi,p / 1 )5(k 1 -k' 1 ) 

and 

W{P, P') = W {0) (P, P') - r 2 W (2) (P, P') + t 4 W (4) (P, P') - .., (14) 

W«(P,P') = fig // (p 1> p' 1 )^ (,) (P,P / ) (15) 

where in the last case we have factored out the small parameter r 2 to 
show explicitly that the last term in (13) is as small as r 2 . 

The structure of the Hamiltonian (13) suggests now the natural 
perturbative scheme for solution of Schrodinger equation (6). Let us 
represent the wave-function iff as power series in r 2 

^ = ^ + T 2 ^ 2 + r 4 ^ 4 + ... (16) 

then inserting (16) into (6) one immediately gets the recursive set of 
equations for i.e. 

(H + V a + V e c ff )* = EV , (17) 
i-i 

(hq + v s + v e c ff ) if 2l = Eif 2l -Y,{-i) l - s w { ~ 2l - 2s ^ 2s , 1 = 1,2, 3, ... 

s=0 

(18) 

The scheme (17, 18) has a remarkable property, namely, the solution of 
the three-body problem in framework of this scheme can be obtained in 
terms of solutions of two-body problems. Indeed, equation (17) allows 
the separation of variables, so that its solution is reduced to the solu- 
tion of respective two-body equations, at the same time the solution 
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of inhomogeneous equations (18) can be obtained in terms of Green's 
function of the operator Hq + V 8 + V£jj which can be constructed from 
two-body Green's functions due to separability of variables. 

For the specific case of three particles in the continuum we are 
considering the solution scheme outlined above yields the following 
results. The solution of the first equation (17) reads 

^o(P,P m ) = ^ r (k 1 )^„( Pl ), (19) 

where initial state momentums kf 1 and p|" are related to the energy 

E by the formula P*" 2 = kf 2 + pf 2 = E. The function ^ c 1„(pi) 

Pi 

is the momentum space Coulomb wave function corresponding to the 
potential v°ff(jpi, p^). The inverse Fourier transform of ip ce in (pi) we 

will denote by ip ce i n (yi) which explicit form can be obtained from (4) 
Pi 

when n = n e ff and replacing k|™, xi by p|™, yi, respectively. The term 

fc„ (ki ) is the scattering solution to the two-body Schrodinger equa- 
k i 

tion with the potential u s (ki, k^) and is conventionally represented as 
the solution of the Lippmann-Schwinger integral equation 

fe(k) = 5(k - k m ) - ^ _ J n2 — - J d q v s (k, q)fe (q). (20) 

The solutions of inhomogeneous equations (18) are given by recursive 
formulas starting from 4?o(P,P tn ) 

l-i 

*2i = -Gs,eff(E + iO) £(-l)'-^ 2i - 2 ^ 2s . (21) 

Here the kernel of the operator G s>e ff(z) = (^Ho + V s + V^f f — z^j 
is represented through two-body Green's functions g s for potential v s 
and g c e jj for potential v c e ^ by the counter integral 

G s , e// (P,P',z) = ^- £dC g s (^iX,CWeff(Pi,Pi, z - C) 

with counter C encircling the cut of g s in anticlockwise direction. 

So that, we have constructed the formal solution to the CF trans- 
formed Schrodinger equation (6) for the system ppe in the continuum. 
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The configuration space wave function which obeys (3) can be ob- 
tained from ^ by CF transform (5) 

*(X,P in ) = /dP* c0 (X,P)*(P,P™). (22) 

Now one can see, that the structure of our solution (16) and respec- 
tive series in r 2 for ^(X, P m ) generated from (22) by (16) and the 
structure of the representation for the Hamiltonian (13) outline the 
framework of our approach as an alternative to Born-Oppenheimer 
one. It is worth mentioning here, that the formalism given above 
is rather general and with minor evident modifications is applicable 
for the three charged particle systems for the case of different masses 
when the mass of one particle is significantly smaller than the masses 
of others. 

Before proceed further, let us give some explicit formulas for ap- 
proximation to the wave function \I/(X,P m ) generated by our com- 
plete formal solution which we will use in the next Section discussing 
some astrophysical reactions. Introducing (19) and (21) into (22) we 
get 

*(X,P m )=V+ r (x 1 )^(yi) + ^ 2 (X,P m )+0(r 4 ) (23) 
where 

^™( Xl ) = / dki^ 1 (xi)^„(ki) 

and ^2(X,P m ) is given by transform (22) of ^2 calculated through 
^0 from (19) by formula 

^2 = -Ge,eff(E + iO)W^^ . (24) 

The complete investigation of properties of the solutions to (17,18) 
is out of the scope of this paper and will be made elsewhere. Below 
in this section, we consider two points which plays the key role for 
the formalism, namely the singular structure of operator W and the 
structure of correction term ^2 (and consequently ^2) which possesses 
the most important properties specific for all correction terms ^21- 

The kernels of the operators in (15) are given in terms of 

singular integrals (10). These integrals obviously can be computed in 
terms of distributions (generalized functions) as it could be seen from 
(10) in the trivial case / = which yields L(°)(P,P') = <5(ki - k'J. 
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The general case of arbitrary I in (10) is considered in the Appendix 
where it is shown that the most singular part of the integral L l (P, P') 
has the form 

LW(P,P') = 5(k,k')5^(k - k')(k,p- p'} l £V\k,k'). (25) 

Here and in what follows we omit subscript 1 from definition of mo- 
mentums and coordinates using, for instance, k instead of ki and so 
ones. Delta-function 8(k, k') on unit sphere S 2 = {k : k = 1} and l-th 
derivative of delta-function 5^ l \k — k') are defined by 



dk<5(k,k') 5 (k , ) = /(k), 

s 2 

/oo 
dk' 8®(k - k')g{k) = (-1)V°W- 
-oo 

The function C^ l \k, k') being a smooth function of k and k' for /-even 
has the form 

2 



; i( CTO (fc)-ao(fe'))+^ (2fc)- i T(2A: , ) i7 'r(l - i( 7 - t 7 ))' 



Above formulas for L^)(P, P') can be used to compute the action 
of operator on 

Wm (P,P in ) = ±h(k,k in )l2(p,p in ) (26) 



where 



and 



/i(k,k m ) = J dk'6(k,k')6^\k-k'^ 2 \k,k')4j+ n ^') 

h(p,P m ) = J dp , ^) e c // (p,p / )(k,p-p , ) 2 ^n(p , )• 

Both integrals Ij are singular distributions. To make its structure 
more transparent let us find its most singular parts. For the first 
integral it means leaving in the integrand the most singular part of 
Y£ in (k') from (20), i.e. <5(k - k in ) which yields 

/i(k,k in ) = 5{k,k in )5 {2) {k - k in )£ {2) {k,k in ). (27) 
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For the second integral, it is useful to make a linear change of variables 
to get 

h( P ,P m ) = ^rj dq(k,q) 2 vg„(q + p) 

and then using Fourier transform for tpp n (p) rewrite this integral as 
/ 2 (p,p m ) = ^rj dyD(y,k)e~^^ m (y). (28) 

Here the function D(y, k) is given by 

D(y,t) = lim o ^ / dqe^)-^(k,q) 2 . 

In the Appendix it is shown that the main singular part of D (y, k) is 
proportional to delta-function, i.e. 

The latter gives for the integral / 2 (p, P m ) 

/ 2 (p,p m ) = iV e// ^(0), (29) 

N = U eff 1 

e// 2vr 2 3(2vr) 3 / 2 ' 

Now, inserting formulas given above in (24) we can represent the 
correction term ^ as the integral 

* 2 (P, P-) = ~ J dP' £ dCUK k', CWeffiP, p',E-( + iO)x 

(30) 

7 1 (k',k-)7 2 (p',p-). 

This general formula can be simplified if we take instead of full Green's 
function g s (k, k',£) its main singular part which is Green's function 
of the two-body kinetic energy operator 5(k — h')(k 2 — C) -1 - This case 
in fact has the particular physical sense, since taking into account 
that we left only delta-function for ipf the resulting approximation is 
exactly equivalent to the neglect of the short-range potential V s from 
the very beginning. Formula (30) is transformed in this case to 

tf 2 (P,P™) = (31) 
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~N eff h(k,k in ) J dp'g c eff (p,p',E - k 2 + iO)^„(0). 

The configuration space representation for ^2 which can be ob- 
tained from formula (22) is reduced now to the integral 

^ 2 (X,P m ) = 

-^e//^n(0) J dk^^)h(Kk m )9eff(y,^E - k 2 + i0). 

Final form for this integral follows immediately from delta-functional 
factors of Ii and reads 

^ 2 (x,p m ) = -iiv e//V; - n(0)x (32) 



^ \t 2 C^ (t, k m W^ n (x)g c eff (y,0,E-t 2 + iO) 
01 t=k" 

This formula describes the correction term ^2 for the ppe system when 
strong pp interaction is neglected and at the same time is approxima- 
tion to the term ^2 in the general case. 

The formula (32) is useful for constructing the coordinate asymp- 
totics of \I/2(X,P m ) as y — ► 00. One needs to use well known coor- 
dinate asymptotics of Coulomb Green's function as y — > 00 and y' is 
bound 

y'. ^ + ») ~ exp( " s, ";| I ' og2si,} ^(y') 

which can be found for instance in [7]. This asymptotic formula gives 
the following asymptotics of ^(X, P m ) 
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^ 2 (X,P iri )~ (33) 



exp \ip m y — i^~r log 2p m y\ ( ( u 

A(x, k«\ p-, y) ^ [X + O [y k - 

where the amplitude A has the explicit form 

„4(x,k"\p"\y) = -±N eff r p in(0)x 



in 



:tn 



8P_ r 

dt 2 



t=k v 



Here by O (y^r^j we have denoted terms corresponding to derivatives 
of exponential factor in (32). The order of this terms shows the range 
of validity of the asymptotics (33), i.e. y^ has to be small, what 
in terms of masses must be equivalent to the fact that yr 2 has to be 
small. Let us emphasize that all treatment above devoted to consider- 
ation of three body charged particles at kinetic energies, comparable 
with corresponding potential energy, which means that Coulomb in- 
teractions are essential. In this situation for the systems, consisting 
from heavy and light particles, one can develop the adiabatic descrip- 
tion, which actually means smallnes of momentums of heavy and light 
particles ratio This smallnes obviously appear due to the small 
mass ratio parameter r, introduced above. 

Therefore in each fixed order of expansion in small parameter r 
one should neglect also by all terms, proportional to the ratio 



3 Astrophysical Examples 

Now let us discuss ways of describing some reactions of the pp-cycle 
on the Sun, which can be done on the ground of 3-body wave function 
given by (23). In other words, we will consider situations when in 
the initial state the system consists of three charged particles in the 
continuum and the mass of one of them is much smaller than other 
masses. 

The first example gives the reaction 

p + p + e ->■ d + v (34) 
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considered in [6] . As it follows from the form of the main term in the 
right hand side of (23), with very good accuracy we have separation 
of the Jacobi coordinates in the wave function of the initial state for 
the reaction (34). This means that the rate of three-body process (34) 
can be expressed in terms of a binary process 

p + p^d + e + +u (35) 

This is just the main result of paper [6]. Now it becomes clear that 
the physical background of the above result from the point of view of 
the few-body theory consists in two points: 

a) the system has two heavy and one light particle such that the 
parameter jfr <C 1 and therefore one can neglect the second term in 
the right-hand side of (23). 

b) heavy particles are slow enough to neglect higher partial waves 
in their relative motion. One should emphasize that free "effective 
charge of the initial nuclear system" Z, introduced in [6], can now be 
fixed to value Z = 2 which is supported by the structure of (23). 

Let us consider another example of 3-body initial state 

p+ 7 Be + e, (36) 
which can produce 8 Be or 7 Li nuclei via the following reactions 

/ 8 -Be + e( 7 ) 

p + 7 Be + e 

\ 7 Li + v + p 

First, from the previous discussion one can see that due to differ- 
ent masses of heavy particles in this case the contribution from the 
linear term r is nonzero in contrast to the p + p + e system, and this 
contribution should be estimated. 

If the electron in state (36) is in the continuum, then again due to 
the separation of Jacobi coordinate in the first term of (23) the rate 
of the proton capture from the initial three-body state (36) can be 
expressed via the rate of the binary reaction p + 7 Be B + 7. 

However, the rate of the electron capture from the initial three- 
body state (36), as it follows from (23), (modified for the state (36)), 
will be defined by the Coulomb wave function of the electron moving 
in the Coulomb field with the charge Z = 5 instead of Z = 4 for the 
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capture from the two-body state 7 Be+e. This means that the produc- 
tion rate of 7 Li from the three-body state (36) cannot be expressed 
via the binary (e+ 7 Be^ 7 Li+^) reaction rate. Roughly speaking, the 
ratio of these rates will be proportional to the ratio of the correspond- 
ing electron Coulomb functions at energy in the center of the Sun E s . 
In other words, 



w 3 

W2 



M0,E S ,Z = 5)\ 2 



V> C (0,£ S ,Z = 4) 



Now let us discuss the problem of screening of the Coulomb in- 
teraction between two protons for the system p + p + e. We restrict 
ourselves by lowest order in the ratio m e /m p for the three-body wave 
function, i.e. consider only first term in the (23). It is evident, that 
the screening effect in this approximation appears due to the elec- 
tronic wave function ip c (p, y), where y = \J mT+2m^ (t + r ) > being 
the distance between protons and r being the distance between elec- 
tron and one of the protons. Taking the asymptotics of this function 
in the region where R 3> r, one can see that the Coulomb phases of 
pp wave function and electronic wave function can cancel each other 
for the specific configurations of initial momentums k and p of three- 
body system under consideration. Hence the resulting motion of two 
protons in this configuration would be described by plane wave, which 
means the total screening effect. 



4 Conclusion 

In conclusion we emphasize, that the CF-transformed three-body Hamil- 
tonian (13) for the system of two heavy and one light particles can be 
used for realization of adiabatic expansion which is alternative to the 
Born-Oppenheimer one. This approach allows to treat screening ef- 
fects on the microscopic level. In astrophysical examples it was shown, 
that in the lowest order of small parameter r, it is possible for some 
reactions only express rates of 3-body processes in terms of binary 
ones. 
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5 Appendix 



In this section we give a brief description of evaluation stages of inte- 
grals LW(P,P') and £>(y,k) which are defined as 

LW(P,P') = A Iim ) LW(P,P',A), 
L«(P,P',A) = J dx ie - A ^Vr i (xi)(x 1 ,p 1 -pi)V^ 1 (xi) (37) 



and 



£>(y,k) = lim £>(y,k,A), 

A — ►H-u 



D(y,k,A) = | dqe- A ^^(k,q) 2 . (38) 

The technical tool we use for calculations of prelimiting integrals (37, 
38) is so called weak asymptotics. For Coulomb wave function ^£(x) 
as y — > oo [7] it reads 

r/- exp (— ifcx — zwo(x, /c)) ,, „ . exp (ikx + iwo(x, k)) 
<5(x,-k) s c (fc,x,k) 



x 



where wo(x, k) = —7 log 2kx, 7 = n/2k and s c (k, x, k) is the Coulomb 
s-matrix 

S c (fc,X,k) = 2Z7T |x-fc|W 

uo(/c) = argT(l +17). If n = this formula is reduced to the weak 
asymptotics of plane wave [9] 

<^ ,M ~ - ^>«} . 

(40) 

The last key formula we need to compute above integrals is 

/<i, c (M,%(k) = .-,(*) + tl^L^^^L, 
with 

/■27T 

G(0)=2^(x), G(t)= / 
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In these equations, 9 and (f> are spherical angles and t = 1 — cos 9. In 
the sense of distribution it means that 



s c (£;,x,k) = e 2wT °<5(x,k) + ... (41) 

The integral from (37) for A / can be calculated by following 
procedure, inserting (39) into the integral, leaving the main singular 
part of s c from (41), and using the formula [10] 



/ 

Jo 



" j- Jlo, ±itx-\x _ r ±W2=FW2 r ( l + la ) 

o (t ± iA) 1+M 



one arrives at 



L«(P,P',A) = -i— <J(k,k , )e- i<J() ( fc ) +CT °( fc ')x 
<-k, P - p') l B t (k, k'X) + (k, p - p') l Bf(k, k'X) 



where 



Bi(k, k'X) = e ^oW-.a,(^)( 2fc) -7 (2fc0 .7 e ^ r(l + I - »(7 - y)) _ 

It remains to evaluate the A — > limit, which can be done by means of 
the following representation which one can verify by straightforward 
calculations 

x^+oj-oo (t ±iA)"+v T(n + i/i)l p y vn J 
/- 1 dt^-^^-g^- 1 )^)] ( r 1 [°°\ dtg<- n -V(t) \ 

y_i (* ± io) 1 ^ + u-oo A / (* ± io) 1 ^ j ' 

Using symbolic notations and leaving explicitly only main singular 
part 

^-tM^-^-ljn-l^n-l)^) [1 _ e ±^] + ... 



(i ± i0) n +^ T(n + i/x) 

we arrive finally at the formula (25) given in the main text. 

By very similar way the integral (38) can be evaluated, i.e. usage of 
(40)which is at the same time the asymptotics of plane wave as q — ► oo 
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helps to calculate the integral over q and then the radial integral over 
q gives 

ly ' ' j (2vr)3/2 (y2 + A 2 )2 • 
The A — > +0 limit is now straightforward 

lim D(y,k,X) = —^—5(y). 

A^+o VJ ' ; 3(2vr) 3 / 2 y ' 
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